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A two-dimensional gas of non-interacting quasiparticles in a nearly periodic potential and a 
perpendicular magnetic field is studied. The potential is a superposition of a periodic potential, 
induced e.g. by a charge density wave or a vortex lattice of a type-II superconductor, and a weak 
random potential due to impurities. Approximating this model by Dirac fermions with random 
mass, random energy, and random vector potential, we evaluate the density of states and the Hall 
conductivity using a self-consistent approximation. We obtain a singular average density of states 
p{E) ~ po + l^-l", where a decreases with the strength of randomness. The Hall conductivity has 
a plateau which is destroyed for strong tunneling through the saddle points of the nearly periodic 
potential. 

PACS numbers: 71.10.Ca, 72.20.-i, 73.40. Hm 



I. INTRODUCTION 



A two-dimensional system of non-interacting quasiparticles is considered in a periodic potential with a strong 
perpendicular magnetic field. The periodic potential can be understood as a pinned charge density wave, formed in a 
two-dimensional low-density electron gas due to Coulomb interaction. This model is related to a possible formation of 
a Wigner crystal in a two-dimensional electron gas with a perpendicular magnetic field, a problem which has attracted 
considerable attention recently |^J^ . Another possible realization of the model is a system of quasiparticles in a layered 
type-II superconductor in the presence of a magnetic field, where a periodic structure due to the vortex lattice in 
the Abrikosov phase arises. Due to impurities, the periodicity of the potential is not perfect in a real system, and we 
have to incorporate weak random fluctuations of the potential. Therefore, the total potential is a superposition of the 
periodic potential and weak randomness, creating a nearly periodic potential. 

We start with a continuum model which we map to the network model, following the ideas of Chalker and Codding- 
ton 1^ for the quantum Hall system. The resulting model is defined by an evolution operator for the quasiparticles 
on a two-dimensional lattice, where the latter is an approximation of the nearly periodic potential. Two types of 
scatterings are possible. One describes the scattering along the equipotential lines, while the other corresponds to the 
scattering through the saddle points of the potential. This model undergoes a transition between Hall plateaux of 
the integer quantum Hall effect, the quantum Hall transition (QHT), when the scattering through the saddle points 
exceeds a critical value. The network model has been very successful as a starting point for numerical studies of the 
QHT [H^J^,^ . Unfortunately, the model is not easily accessible to analytic methods due to its internal lattice struc- 
ture. However, a large scale approximation, for weak disorder, sufficient for the investigation of transport properties, 
has been applied to the network model by Ho and Chalker ||^ . They found that the network model is equivalent to a 
model of Dirac fermions. This result indicates that the network model is similar to the tight-binding model with half 
flux quantum per plaquette, where the large scale approximation also leads to Dirac fermions ||,^,|l^. However, the 
corresponding randomness in the resulting Dirac fermions is different: the tight-binding model with random potential 
leads to Dirac fermions with a random mass whereas the network model of the nearly periodic potential leads to Dirac 
fermions with a random vector potential and a random energy. 

Dirac fermions with randomness can be treated by field theoretical methods J0,|ll|,|l2| . A particular case is the 
purely random vector potential, which can be solved exactly by bosonization P,p^. Other types of randomness, to 
the best of our knowledge, cannot be treated in this way. Therefore, we will apply a self-consistent approximation 
scheme. In order to control our approximate calculation we will compare it with the exact result of Refs. P,p^. 

The paper is organized as follows. In Sect. II we introduce the continuum model, and derive a network model of 
quasiparticle loops, which is equivalent to the network model of Chalker and Coddington. In Sect. II. A we briefly 
recall a tight-binding approximation of the periodic potential [9|Jl^ . The large scale approximation of the network by 
Dirac fermions is discussed in Sect. II. B. As physical quantities, we define the average DOS and the Hall conductivity 
axy in terms of the Green's function. The main part of our calculation is presented in Sect. HI, where a self-consistent 
approximation for the average Green's function is worked out, as well as its consequences for the average DOS and 
the Hall conductivity. Finally, we discuss our results in Sect. IV. 
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II. THE MODEL AND RELATED PHYSICAL QUANTITIES 



We consider the dynamics of a particle with charge e in the potential V{x,y) — Vb cos(7raa;) cos(7ra?/) + SV{x,y), 
where dV{x,y) represents a weakly fluctuating disorder potential. We start with a discussion of the classical motion, 
taking into account the Lorentz force created by the perpendicular field B = (0, 0, B): 

where r = {x,y,0) is the coordinate vector. For a sufficiently small mass M the first term can be neglected, and the 
resulting first order differential equations show that the particle follows the equipotential lines of the potential: 



d_ f x\ ^_^{ dV/dy 
dt\yJ^eB\ -dV/dx 



(2) 



The approximation by the first order differential equation is equivalent with neglecting the cyclotron motion of the 
particle and reducing the dynamics to the motion of the guiding center. For _B > a particle with energy E ^ Q 
follows closed counterclockwise orbits around a hill {E > 0) or clockwise orbits around a well {E < 0) (Fig. 1). 
Only a particle with energy E ^ can travel through the whole two-dimensional space because there are extended 
equipotential lines at a; = {2k + l)/2a or at y = {2k + l)/2a {k = 0, 1, ...). The behavior near E = can be quantized 
using a lattice approximation. This approximation is based on the orbits of quasiparticles at a given energy E which 
form a network of loops. Quantum effects are described by tunneling between these loops in regions where the loops 
are close to each other, i.e. at the saddle points of the potential. For E < the tunneling occurs between the loops 
in Fig. 2 whose centers (minima of the potential) are located at {x,y) and {x ± l,y ± 1). For i? > we have a 
complementary situation in which the preferred loops appear with the opposite current direction. Apart from this 
difference, the subsequent analysis can be applied to this case as well. We will consider E < in the following. 



A. Lattice Approximations 



The tunneling can be expressed by discretized scattering processes of effective quasiparticle positions on the centers 
of the loop edges, enumerated by 1,...,4 in Fig. 2. The scattering is characterized by a scattering parameter /3 
(0 < /3 < 7r/4), determined by Vb and E. Regarding the loop at {x,y) in Fig. 2, there is a quasiparticle state 
^{t) = (4'i(i, X, y), 5'3(t, X, y), 4'2(i, x, y), 5'4(i, x, y)) at time t whose components are associated with the edges j on 
the loop. At the later time t+1 this state can be scattered along the loop, or can tunnel through one of the four saddle 
points to the loops at {x ±l,y ± 1). For example, the component ^'i(i, x, y) can be scattered either to ^2(^ + 1, a;, y) 
with rate C = cos/3, obtaining a phase (j)2, or can tunnel to ^'4(t + 1, x + 1, j/ + 1) with rate S — sin/3, obtaining a 
phase 04. All possible processes can be combined as follows: 

+ (n ^)*(*)' 

where 

and (ttl.), for example, are lattice shift operators in the positive (negative) x (y) direction. The phases 04 
are associated with the centers of the loop edges. The weak random potential 6V induces weakly fluctuating phases 
with = for j = 1,2,3 and (^4) = tt. This choice ensures the presence of half a flux quantum per current 
loop on average. (Notice that the original network model |^ has strong disorder due to a uniform distribution of the 
phases on the interval [0, 27r). It is not clear if our model with weak disorder belongs to the same universality class. 
Since strong fluctuations of the phase correspond to strong fluctuations of the external magnetic field, the fact that 
a random magnetic field is a relevant perturbation indicates that both models may be qualitatively different.) 
The random potential SV also affects the tunneling between the loops. This induces randomness in the lattice shift 
operators t^'^. The latter are statistically independent for different nearest neighbor pairs on the lattice. We believe 
that /3 > 7r/4 is not realistic for the physical model because this describes a situation where the quasiparticles prefer 
to tunnel through the saddle points rather than go along the equipotential lines. 
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The (discrete) time evolution of the quasiparticle state on a loop can also be described by the evolution operator 
W as *(t + 2) = W*(i) §, with 

fo mV {MN \ ... 

Assuming (discrete) translational invariance (i.e. SV = in the original model), we can diagonalize the scattering 
matrix by applying a Fourier transformation (x, y) — > (fci, fc2). Thus we find the eigenvalues of W, 

Ai/2(fci,fe) = 2S'Ccos fci cosA:2 ± Wl - 'iS^C'^ cos^ h cos^ fca. (6) 

The corresponding dispersion relation, i?x/2(fci7 ^2); obtained via X1/2 = exp(— ii?]^/2)j has a gap around E = 
depending on /3 (see also 0). The gap vanishes at fci = ^2 = 0, ±7r for (3 = tt/4, i.e. for C = 5 = 1/2. Thus the 
quantum model reflects the behavior of the classical model, where infinite equipotential lines exist only at zero energy. 
The effect of disorder on this behavior will be studied subsequently. 

We briefly discuss an alternative lattice approximation of the almost periodic potential. It is based on the tight- 
binding representation of the continuous model. For this purpose we regard the minima, V{x,y) = —Vq, and the 
maxima, V{x,y) = Vq, of the periodic potential as a lattice. Quasiparticles can hop between neighboring lattice 
points with hopping rate t, i.e. for hopping between a maximum and a minimum. Moreover, it is reasonable to 
take the hopping between next-nearest neighbors into account also which are either pairs of maxima or pairs of 
minima, and associate a ho ppi ng rate t' with this process. The magnetic field leads again to a Peierls phase factor in 
each of the hopping terms Q. Disorder can be introduced by adding random fluctuations SV to the lattice points. 
This lattice approximation differs from the above network model because there are two types of hopping now. The 
nearest neighbor hopping between a maximum and a minimum has been excluded in the network approximation, 
since tunneling is considered only between states at the same energy. 

B. Large Scale Approximation: Dirac Fermions 

We are interested in the properties of the quasiparticles on large scales. Therefore, we consider the large scale 
approximation of the evolution operator W by expanding the non-local part in a Taylor expansion t^/(r) = i^(r)/(r± 
aex) ~ t^(r)[/(r) ± aVa;/(r)]; for simplicity we set a = 1 in the following. The evolution operator then reads for 
P = 7r/4-F5/3, with \Sf3\ < 1, 

- 1 + ( -2%+vV-\ ) + 

and 

- 1 + ( 2^ - vi + & 11 - vr;tif ) + --d-- '^""^ 

with the effective random vector potential components — (^i — <j)^)/2^ Ay — ((f>4 — tt — (j)2)/2, and with — 
(</>! -I- • ■ • + 04 — 7r)/2. We can express the results (^, (|^) using two-dimensional Dirac Hamiltonians Hj^i^ and Hnm, 
defined through MN « 1 — iHmn and NM w 1 — iH^M, where 

Hmn = -2(5/3(72 - («Va; + Ax)a3 + (iVy + Ay)ai + Hmn 
Hnm = 25(3(72 + {iVy + Ay) + A^)(Ji + Hnm (9) 

with Pauli matrices aj. The terms Hmn, H^m are contributions from the random operators t^'^. Here m = 25(3 
appears as a Dirac mass. Since 5(3 < 0, the Dirac mass cannot be positive. Moreover, it is convenient to rotate the 
2x2 matrices Hmn and Hnm such that the mass term is in the diagonal. This leads to 

Hmn = - (iVa; + Ax)ai - {iVy + Ay)a2 + Hmn 

Hnm = maa + {iV^ + A^)ai + {iVy + Ay)a2 + Hnm- (10) 

The random terms, related to the random phases <j)j and the random factors t^'^, are expanded in terms of Pauli 
matrices as 'Y^i=o'^i^i, where I — \^2 are contributions to the random vector potential, Z = is the contribution to 
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the random energy, and I = 3 the contribution to the random mass, (ctq is the 2x2 unit matrix.) The corresponding 
expansion for Hj^m is given by X];=o '^i^i ■ assumed that Vi and V/ have a Gaussian distribution with zero mean 
and correlations 

{Vl,rVl>y)v = {Vi[rVl'y)v' = Srr'Su'9l, (11) 

where r,r' denotes the lattice points. In principle, Vi and V/ are correlated. However, these correlations do not play 
a role here because the Green's function is block-diagonal with respect to Vi and V/ . 

We emphasize that a non-compact continuous symmetry exists in the case — — (pure random vector 
potential) and m = 0, namely 

HmN = [(1 + C')'/'^To + C^3]HMN[il + C')'/''To + C^g] (12) 

and similarly for Hmm, with — cxo < C < oo. This symmetry indicates that the pure random vector potential is 
qualitativly different from the general case with go + 33 > 0, e.g. the symmetry is reduced to a discrete one with C = * 
for averages of physical quantities li m — E = Q. 

The time evolution of the quasiparticles in the network is given by iterating the time step operator W. This requires 
the evaluation of (W")\/ for large values of n, since we are interested in the long time behavior. We Fourier transform 
W" first, and later average with respect to the random potential V . In a first step, 

^ W"e""' = (1 - We*'^)-\ (13) 

where the frequency lj has an infinitesimal imaginary part e > 0: lo — E + Now we return to the Dirac Hamiltonian 
to determine the large scale behavior. For small frequencies we obtain the Green's function of Dirac fermions 

(l-We-)-i« f*^*'^^"''' , \ \^G{u:). (14) 

The evaluation of {G{uj))y in a self-consistent approach is described below (Sect. III). 



C. Density of States and Hall Conductivity 

According to standard Green's function theory, the average DOS at a given energy E is obtained from the Green's 
function G as 

p{m, E) = - lim(Im[G'ii rr{E + ie) + G22 rr{E + ie)])v- (15) 
The DOS of a pure system [gi — 0) is linear around = ^ and vanishes at _E = for all values of m. 

The current density in a Dirac model can be calculated directly from the Kubo formula |^ or from the response to 
an external static vector potential qy [Q. The effect of qy is a change of the boundary conditions in the y-direction, 
a concept extensively used in numerical investigations of Anderson localization to study the existence of delocalized 
states [|l^. The response to the vector potential gives the Hall conductivity a.j:y in terms of Green's functions |p|,p^ 

O-^y = / ^ (Tr2|crj;(iJMAf - iz - E)^1,{Hmn - iz - E - qyCry)^^A 

+Tt2{(jAHnm -iz- E)-1,{Hnm ^iz-E- qyCTy)-,]^Y) (1^) 
This gives, in the absence of disorder and with periodic boundary conditions, and for qy 0, this leads to 

a„j = ^sign(~m)e(|m| - |£;|), (17) 

where O is the step function. Thus the Hall conductivity has the plateau value /h at £' = 0, since m < for 
(3 < 7r/4, while it is undefined at m = 0. This behavior of the transport quantity axy is unphysical, and will be 
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replaced by a continuous behavior when we take disorder into account: then the DOS is found to vanish at -E = for 
all m, which indicates that the Hall current flows on the boundary. Therefore, the result ( |l7|) depends strongly on 
the choice of the boundary conditions (|J^; conditions other than periodic, can suppress axy For the tight-binding 
model of Sect. II. A, axy was determined for E ^ ||l^ to be 

axy^^eif'^V,'). (18) 

Comparing this with the result (|l^), we see that the network approximation is similar to the tight-binding approxi- 
mation provided t'^ > Vq. 

For the special case of a random vector potential we can evaluate the dissipative conductivity and the localization 
length in a simple way. Using the identity (which holds for m = 0) 

Gi-uj) = -a3G{io)a3, (19) 

we may express the two-particle Green's function {G{u))G{u}*)) by the product of one-particle Green's functions on 
the same complex half-plane. As a consequence of this property, the critical properties of the two-particle Green's 
function are identical with those of the one-particle Green's function. This explains the finding that the critical 
exponent of the localization length is identical with the decay exponent of {G{uj)) in an explicit calculation [Bj. 



III. SELF-CONSISTENT APPROXIMATION 



The full Green's function {G)v can be related to the Green's function Go in the absence of randomness, 

Hmn)v 
Lo - {Hnm)v 



Go(..)=("-^^*^^^^ V ) , (20) 



through Dyson's equation: 



{G)v = Go + Go^{G)v. (21) 



The self-energy S, in lowest order perturbation theory, is given by E ~ {VGqV)v ■ Expanding S in terms of Pauli 
matrices, S = Eo + Sio"i + ^2(^2 + Sscts and, taking into account (G)y^ = Gg ^ — S, we see that Ei and E2 shift 
the gradient operators in {Hmn)v and {Hnm)v- Since the parity is conserved for Go as well as for {G)v, these 
contributions must vanish. The remaining contributions, So = — «'7 and S3 = —mg, can be evaluated in a self- 
consistent approximation, where we use second order in perturbation theory, and replace Go in the expression for S 
by (Gg ^ - S)-\ with the result 

E=(y(Go-i-I])-V)y. (22) 
With gij = gi + gj, this self-consistent equation leads to two coupled equations, 

(.912 - 903)1 

ms = m- — (23) 

1 - (912 - 903)^ 

and 

7? = (77 - iuj)gl (24) 



with the integral 



tPk 

[(m + m,y + (r; - lu^f + k^]-^—— 



1 / \ 

= — In 1 + . (25) 

ZTT V (m + ras)" + [rj — luy / 

We introduced A as an ultra-violet cut-off, and g = go + ■ ■ ■ + 93- 

For w = 0, there are two solutions of (p^: 77 = 0, and a second one with r] ^ and / = g^^. The latter implies 
with (psh a renormalization of the Dirac mass. 
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m = m + Trig 



25, 



'03 



and with (Ea), we find 



(26) 



(27) 



The multiplicative renormalization of the Dirac mass in Eq. ( p6| ) diverges as go3 ~^ 0. This indicates a special behavior 
for a pure random vector potential. Although this case is not interesting in terms of the network model, it has been 
discussed extensively in the literature P jTTlJli ] . 

For TO = (oj 7^ 0) we obtain with Eqs. (|24), ( |25| ) an equation which determines 77, 



Given rj, we can evaluate the DOS as 



(ji — iujy 



(28) 



p(m = 0,E) = -Re(?7), 
9 



which for E = vanishes if > to^, where 



nir — ± 
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(29) 



(30) 



where we set A = 1 for simplicity. In this case the DOS as a function of to follows a semi-circle (Eq. (|2^)), multiplied 
by (7/2^03, with the maximum at to = given by p(m = 0, -B = 0) = e^^^^/g. The radius of the semi-circle is 
2(503/5)6"'^^^: and thus vanishes for go3 ~* 0- This means that the DOS vanishes for a pure random vector potential 
(503 = 0) for any to 7^ 0, but jumps to p{m = 0) = e~^^^^^ / gi2 at m = 0. This again indicates a singular behavior 
(instability) for the pure random vector potential, at least in terms of our self-consistent approximation . Such 
behavior was also found in perturbation theory around the (integrable) case (random vector potential only), where 
infinitely many relevant operators appear due to the contribution of the random Dirac mass pT| . 

The density of states p{m — 0, E) behaves smoothly for E Q, a.s shown in Fig. 5. It shows an effective power law, 
e~^/9 _|_ ^ith a cut-off at large E (depending on A), except for very small and very large energies. As can be 
seen in Fig. 5, the exponent decreases, starting at a = 1 for g = 0, with increasing g. 



The Hall conductivity can also be calculated within the self-consistent approximation. For this purpose we replace 
the Green's function in Eq. (|l|) by (Gg ^ - 1])"^ This gives an expression identical to that found for the random 
Dirac mass ||T^, except that rric is replaced by the new critical mass (30): 



— sign(— to) 1 — (2/7r)arctan(-\/ rri^/m? 



(31) 



This describes a Hall plateau with a continuous decrease to zero near the critical scattering parameter (3c = 7r/4, 
which implies that disorder changes the (unphysical) discontinuity of the pure system (O) to a more realistic behavior. 



However, for go3 = 0; 5 > 0, the critical mass to^ vanishes. Thus (Jj-y is undefined for to = TOj,, similar to the pure case. 
Consequently, the random vector potential alone is not sufhcient to create a continuous behavior of Uxy However, 
since there exist exact calculations for the random vector potential Q and numerical results [T^ , we wish to compare 
our self-consistent approach with those. In order to avoid problems with singularities at to = we assume to > and 
set TO — > only at the end. In particular, we have to choose that solution of Eq. (|4|) which vanishes with vanishing 
frequency lj. Going back to the self-consistent equation (p^), we can express / of Eq. (^5|) by the renormalized mass 
TO as (1 — l/x)/gi2, where x = ifi/m. Consequently, ry - compare with Eq. ( p^ ) - is given by 77 = icij(l — x). From 
Eqs. (H), (H), we obtain 



X = < I- 



I 2tt 



In 



1 



1 



x'^{m? — w^) 



This implies, with lo 



2 < 



|l + |^ln \x\m^-oj^)\] 



(32) 



(33) 
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Assuming weak disorder {gi2 ^ 1) we exponentiate the right-hand-side to get 

X ~ _^2)pgi2/2^ (34) 

and 

x^im^ - uj^)-9i2/[2n(i+g,2M]_ (35) 

The decay length ^ of the average Green's function reads in terms of the self-consistent approximation (in units of 
the lattice constant) 

e= [7fi2 + (r;-ic^)2]-i/2_ (36) 



Thus, from (35) and inserting m, rj ^ iLu{l — x), we find 

e^(m2-c^2)-i/2(i+<„.A)^ (37) 

This agrees with the earlier results |9|,[T3t, obtained from a renormalization group argument and a bosonization 
approach, respectively. In the latter it is not only restricted to weak disorder, but holds true for all values of 512. 

IV. DISCUSSION AND CONCLUSIONS 

Starting from a two-dimensional system of non-interacting quasiparticles in an almost periodic potential with 
magnetic field, we have derived an effective model for large scale properties. The latter is a model of 2D Dirac 
fermions with random mass, random energy and random vector potential. Our derivation is analogous to that of Ho 
and Chalker for a random potential This model shows a gap due to the periodic potentia,l which disappears if the 
rates for clockwise and counterclockwise scattering of quasiparticles are equal. The gap is proportional to the Dirac 
mass. The average Green's function of the random Dirac model has been evaluated in a self-consistent approximation. 
The latter is equivalent to the saddle point (or large N) approximation of the model ||2^ . It gives a self-energy which 
consists of a multiplicative renormalization of the average Dirac mass, m —>■ mg/2gQ3, and a spontaneous creation of 
a complex self-energy. 

In terms of the self-consistent approximation only the combinations go3 a-nd ,912 enter the results. This supports the 
conjecture 0| that the randomness of the scattering rates at the saddle points is irrelevant (i.e. there is no qualitative 
difference between gs > and gs = 0) as long as (70 > 0. In other words, the randomness of the phases is sufficient 
for the qualitative description of the QHT, at least for the DOS and Hall conductivity. 

The mass renormalization is only a factor 1/2 in the absence of a random vector potential (512 — 0), but it grows 
with gi2, and 77 is a function of the renormalized mass and the (real) energy E. The average DOS is proportional to 
the real part of 77. It shows a semi-circular behavior with width nic, given in Eq. (^o|). This width vanishes if only 
a random vector potential is present, indicating a qualitatively different situation. As a function of E, it describes 
an effective power law, unless E is very small. The exponent of the power law a decreases with disorder strength, 
starting at a = 1 for the pure system. In the case of very weak disorder (g ~ 0), we take the g*'' power of Eq. ([2^), 
and expand the result in powers of g, with the result 

V=^§-M[l^XV{E + ^ef]+0{g'). (38) 



This yields with Eq. (29) the well-known linear behavior of the DOS for E < A , 

p(m = 0,£;) = Me(A2-i?2) + O(5). (39) 

The overall behavior of the DOS is in very good agreement with a numerical result for a similar Dirac model p8| . 
Furthermore, a non-zero DOS near the QHT, where axy deviates from the Hall plateau value e^//i, as found by us, also 
agrees with numerical observations |2^,|2^. The Hall conductivity a^y has a plateau in the regime where scattering 
is dominated by states localized on the loops. However, in the absence of disorder, axy is not defined if scattering is 
equally probable along the loops and through the saddle points (/3 = 7r/4). This problem disappears in the presence 
of disorder, since a^y then changes continuously from the plateau value a^y = e^/h to a^y = at l3 — n/A. 

We conclude that the special properties of the pure system are cured by the randomness: the vanishing average 
DOS at = is elevated to a non-zero value. Although this value is exponentially small, it indicates that a band 
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of states exist where the quasiparticle can tunnel with high probabiHty through the saddle points of the periodic 
potential (i.e. near j3 = 7r/4). The tunneling, which can be understood as quantum percolation, obviously destroys 
the plateau of the Hall conductivity At least for the case of a random Dirac mass, it is known that part of this 
band consists of delocalized states leading to a non-zero conductivity Oxx ~ (e^/^7r)/(l + g/'^.i^) Numerical 
calculations indicate that delocalized states also exist in the network model. Moreover, for the random vector 

potential we found a divergent localization length when m? —>■ uP' within our self-consistent approach. All these 
results are strong hints that there are delocalized states for the fully random Dirac Hamiltonian, in the regime of 
strong tunneling through the saddle points of the nearly periodic potential. 
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FIG. 1. Equipotential lines of V{x,y) = Vb cos(7ra:) cos(7r2/) (arbitrary units). 












FIG. 2. Network of equipotential lines for two fixed energies ±E close to zero. The arrows indicate the directions of the 
quasiparticle motion due to the magnetic field B > 0. 



FIG. 3. Tight-binding approximation of the periodic potential, where t is the nearest neighbor and t' the next-nearest 
neighbor hopping rate. 
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FIG. 4. Density of states for weak {g = 0.2) and intermediate {g = 0.7) disorder with scattering parameter f3 = 7r/4. The 
minimum at S = is e~'^^^ /g (see text). 
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